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g ■ ABSTRACT 

I We report the successful application of a simple recursive proportional- 

jy-^ ' feedback (RPF) algorithm to control chemical chaos observed during the 

■ electrodissolution of a rotating copper disk in a sodium acetate/acetic 

acid buffer. The RPF method is generally applicable to the control of 
chaotic systems that are highly dissipative and adequately described by 
^ , one-dimensional maps of a single variable. The successive minima of the 

' measured anodic current generate a return map that is used to implement 

1^ ■ the control strategy and the anodic potential is the control parameter. Un- 

^ . stable periodic orbits of period one and period two are stabilized by applying 

• i-H , small perturbations (~ 0.05%) to the anodic potential on each cycle of the 

^ I periodic oscillation. Experimental evidence is presented to indicate why the 

' RPF method was necessary in this system and the theoretical robustness of 

the algorithm is discussed. 



1. Introduction 

Two experiments were recently reported describing the control of autonomous 
oscillatory chemical systems in the chaotic regime. Petrov, Caspar, Masere, and 
Showalter0 stabilized periodic orbits in the Belousov-Zhabotinsky reaction and we 

t Published in "Proceedings of the 2nd Conference on EXPERIMENTAL CHAOS" edited by W. 
Ditto, L. Pecora, S. Vohra, M. Shlesinger, and M. Spano (World Scientific, River Ridge NJ, 1995) 
pp 304-316. 
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recently reporteda the first experiment to control chaotic oscillations in an electro- 
chemical cell. In controlling the electrochemical cell we found it necessary, in general, 
to use the recursive proportional-feedback (RPF) algorithm recently developed by 
Rollins, Parmananda, and Sherardi. The RPF algorithm is simple to apply and is 
generally applicable to highly dissipative systems that are well described by a one- 
dimensional Poincare return map of a single measured variable. In this paper we 
report the further use of the RPF algorithm to stabilize both period- 1 and period-2 
orbits in the electrochemical system. We also present experimental evidence that 
suggests why we found it necessary to use the new RPF method instead of the simple 
proportional feedback method suggested by Petrov, Peng, and Showaltera and by 
Hunti that often works for systems well described by a one-dimensional return map. 
Finally, we briefly discuss the robustness of the RPF method. This discussion is based 
on an analysis of the RPF method in terms more akin to standard control theory! 
as described recently for the general high dimensional system by Romeiras, Grebogi, 
Ott, and Dayawansa@. This approach treats the experimental system (described by 
a one-dimensional return map of a single variable) together with a linear recursive 
control stratagy applied to a single control parameter as a discrete two-dimensional 
system. The desired periodic orbit is a fixed point of the two-dimensional map. The 
range of acceptable control conditions is determined by the requirement that the fixed 
point remain stable. The RPF control stratagyi makes the fixed point superstable. 

2. Experimental System 

The experimental system was a three-electrode electrochemical cell set up to 
study the potentiostatic electrodissolution of copper in an acetate buffer electrolyte 
solution of sodium acetate and acetic acid. The experimental setup and the gen- 
eral electrochemical behavior of the system is described in some detail by Dewald, 
Parmananda, and RoUinsi'i. Figure |l] is a schematic diagram of the system. 

The anode was a rotating copper disk (5 mm diameter) shrouded in a cylinder 
of Teflon. The supporting electrolyte solution was a mixture of 60 parts glacial 
acetic acid to 35 parts 2M sodium acetate. The potential of the anode was measured 
relative to a saturated calomel reference electrode (SCE) and the cathode was a 
2.5 cm^ platinum foil disk. The emf of the circuit was continuously adjusted by a 
potentiostat to maintain the desired set value of the anodic potential (the potential 
between the anode and the reference electrode). The anodic potential was used as our 
control parameter. Time series data was collected by measuring the anodic current at 
20 ms intervals using a 12 bit A/D converter interfaced to a computer. The period for 
the current oscillations in the chaotic states used in these experiments ranged from 
2 to 4 sec so that the time between data points in the time series is small compared 
to the typical period of the oscillations. The set point for the anodic potential was 
controlled by the computer through a D/A converter with the smallest increment of 
the control voltage being 0.1 mV. 
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Figure 1: Schematic representation of the three-electrode electrochemical 
cell. V is the anodic potential and / is the anodic current. The potentiostat 
adjusts the emf to hold V at the desired set value. 

3. Control Method 

The RPF control strategy!! was used to stabilize both period-1 and period-2 
orbits within chaotic attractors measured by a time series of the anodic current. 
A Poincare section was taken at the time when the anodic current goes through 
a minimum. A one- dimensional return map was constructed from the sequence of 
current minima where In is the current minimum at the beginning of the nth Poincare 
cycle. The value of the anodic potential during the nth cycle is = + SVn. 
According to the RPF algorithm, control is established by adding an increment to 
the anodic potential during the nth cycle given byi'i 



where Ip is the unstable fixed point of the target orbit for the return map obtained 
for a time series taken at ^ = Vq. The value of Ip and the constants K and R, are 
determined from a precontrol experimental procedure carried out using an interactive 
computer program for data acquisition and display. 

The precontrol procedure is described in detail in referencesili. Firstly, the 
fixed point Ip and the slope, /z, of the one-dimensional return map was determined 
by a least square fitting of the displayed return map (/„ versus /„_i) for a sequence 
of current minima collected in the neighborhood of the desired fixed point with the 
anodic potential held constant at = Vq. Secondly, a sequence of current minima 



6Vn = K{In-lF)+RSVr, 



n—li 



(1) 
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were collected and displayed as a return map while the anodic potential is changed 
back and forth between two values; Vn = Vq for n odd and = Vq + AV for n even. 
The value of AV > was chosen to be about the size of the maximum 6Vn used later 
during control. If the system is sufficiently dissipative, as our electrochemical system 
was, the resulting return map consists of two curves called the up and back maps 
formed from alternate (/n-i, In) pairs with n odd and even respectively. The shift 
in the fixed point per unit AV is measured by least squares fitting for the up and 
back maps giving and o;, respectively. The values of K and R were then calculated 
using the RPF relations^ 

K = ^ , R = (2) 

Typically these control constants, K and R, could be determined within 10-20 min- 
utes (corresponding to approximately 200 cycles of the experimental return map) 
from the start of data acquisition. 

With these control constants available, the control algorithm was initiated. The 
anodic potential was changed according to the RPF algorithm of Eq. || whenever a 
minimum of the anodic current, J„, came within about 0.2 mA of the measured fixed 
point Ip. 

4. Results and Discussion 

4.I. Control of period- 1 and period- 2 orbits 

Examples of the successful stabilizing of period-1 and period-2 orbits are shown 
in Fig. H. The rotation rate for the copper disk anode was u = 2670 rpm, the anodic 
potential with control off was Vq = 0.724 V, and the values of the RPF proportionality 
constants were K = —5.0 mV/mA and R = —0.3 for the period-1 control shown in 
Fig. For the period-2 control shown in Fig. u = 2160 rpm, Vq = 0.745 V, 
K = —5.0 mV/mA, and i? =0.21. As described in reference |^, we chose a region 
of parameter space where the electrochemical system exhibits a sequence of periodic 
mixed mode oscillations separated by bands of chaotic behavior as a function of 
anodic potential. The period-1 control was done while in the chaotic region between 
a period-1 state with large amplitude oscillations only and a period-2 state with one 
large amplitude oscillation followed by one small oscillation. We were unsuccessful 
in controlling on the period-2 oscillation in this chaotic band. Successful control of a 
period-2 state was attained in the chaotic band between the period-2 (one large-one 
small) and period-3 (one large and two small) mixed mode oscillations. We found 
that control of the period-2 orbit was much more difficult than the period-1. This is 
partially because the period is about twice as long (about 5 sec) and the feedback 
corrections are made just once each cycle. Thus the system can drift further away 
from its fixed point before the next feedback correction is calculated and applied. 
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Figure 2: The minima in anodic current plotted over a time period during 
which the control algorithm is switched on and off twice. The perturbations 
added to the anodic potential to maintain control are shown in the bottom 
graph, (a) shows period-1 control and (b) shows period-2 control. 



The corresponding return maps for the two cases are shown in Fig. ^ while Fig. ^ 
and ^ show the chaotic attractors and controlled periodic orbits reconstructed from 
the time series of the anodic current using a two-dimensional time delay embedding. 
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Figure 3: (a) The first iterate return map (open circles) obtained from the 
sequence of current minima for the chaotic state shown in Fig. 2a. The 

superimposed filled circles are the minima in the anodic current while the 
control algorithm was implemented, (b) the corresponding second-iterate 
return map for the period-2 case shown in Fig. 2b. 
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Figure 4: (a) The two-dimensional time delay embedding with r = 120 
msec showing the chaotic attractor for the situation shown in Figs. 2a and 
3a. (b) the corresponding period-1 trajectory with control on. 
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Figure 5: (a) The two-dimensional time delay embedding with r = 120 
msec showing the chaotic attractor for the situation shown in Figs. 2b and 
3b. (b) the corresponding period-2 trajectory with control on. 
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4-2. Why use RPF method 



Feedback methods used to control chaos as discussed in this paper are most 
closely related to the algorithm first proposed by Ott, Grebogi, and Yorke (OGY)0. 
The original OGY method required the measurement of two system variables at the 
Poincare surface (in a three dimensional system) to determine the appropriate feed- 
back for control. However, many experimental situations (such as that described in 
this paper) donot easily allow the measurement of more than one variable. It was 
recently foundi0lli that the original OGY method must be modified to apply to the 
situation where single time series data is used. In general, this requires nontrivial 
modification of the algorithm and makes the prescribed change in the control param- 
eter on the nth cycle depend on the changes that were made on previous cycles. It 
was shown in reference ^ that in highly dissipative systems this reduces to a simple 
recursive algorithm where the change on the nth cycle depends only on the change 
made on the (n — l)th cycle as indicated in Eq. 0. Furthermore, it was shown that 
the recursive term goes to zero [R = 0) if the attractor (in the neighborhood of the 
fixed point at the Poincare section) does not shift in the direction normal to its plane 
in state space when small changes are made in the control parameter. 
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Figure 6: Poincare sections reconstructed using time delay embedding of the 
measured time series of the anodic current. The anodic current is minimum 
at times tn and r = 600 msec. The open and closed circles are for values of 
the anodic potential that differ by 4 mV. 
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Figure |^ shows the Poincare section of the reconstructed attractor in a delay 
coordinate embedding for the electrochemical cell at parameter values where we used 
RPF to control on a period-1 orbit. The time tn is the time when the anodic current 
goes through the n-th minimum near the period-1 orbit. Two sets of data are shown 
with the anodic potential, V, held constant in each case. The open circles are for V = 
0.720 V and the closed circles for V = 0.724 V. This figure showing experimental data 
should be compared with Fig. lb of reference While the observed experimental 
shift in the Poincare section is small, there is clear evidence that a shift is present. 



4.3. Theoretical robustness of RPF method 



We have found in our experimental work that there is often a rather large range 
in the values of the RPF proportionality constants K and R that will still be effective 
in controlling the system. This is fortunate since it allows us to make estimates of K 
and R from just a few experimental points in the neighborhood of the desired fixed 
point during the precontrol phase of the experiment. Here, we briefly derive the RPF 
equations for K and R using methods of control theoryQ'Q in order to to shed light 
on the theoretical robustness of the RPF method. 

We start with the general ID return-map equation (Eq. 2 from reference ^ 

Xn+l = fiXn;Pn-l,Pn), (3) 

where the notation is that used in reference ^ X„ is the value of the measured variable 
at the start of the nth Poincare cycle and p„ is the value of the control parameter 
during the nth cycle. For control on a period-1 orbit, we are interested in the natural 
dynamics of the system described by Eq. ^ near the fixed point, Xp, of the period-1 
orbit for p = po; Xp = f{XF,Po,Po)- To first order in 6Xn = {Xn — Xp), Spn-i, and 
6pn, we have 

SX„+l = flSXn + wSpn-l + vSpn, (4) 

where /i = df{XF,Po,Po)/dXn, w = df{XF,Po,Po)/dpn-i, and v = df/dpn- 

In reference 0, the RPF algorithm is obtained by using the optimum possible 
strategy of choosing Spn such that the system is brought to the fixed point as quickly 
as possible, namely two Poincare cycles. Taking SXn+2 = and Spn+i = 0, then the 
first and second iterate of Eq. (^) gives the recursive control algorithm Eq. |l| with R 
and K given by 

K = —- — -, and R = —- — -. (5) 

[flV + W) [jJ,V + w) 

As shown in reference H, = (1 — ^)gb and f = (1 — ii)gui and Eqs. are equivalent 
to Eqs. |. 

Here we take a different approach. We assume there is a bilinear relationship 
between 5pn and 5pn-i) 

6pn = H6Xn + G6pn-u (6) 
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where H and G are constants yet to be determined. Equations ^ and ^ form a 
two-dimensional discrete map describing the full control system including both the 
dynamics of the nonlinear system and the recursive feedback control strategy. The 
two-dimensional map is put in more conventional form if we shift the n labels on 
the pn by 1. This can be accomplished by defining a new parameter, 5p„ = 5pn-i- 
Substitution into Eqs. || and || gives 



5Xn+l = li 5Xn + W6pn + V 5pn+l 
5pn+l = H 5Xn + G5pn- 



(7) 
(8) 



The conventional 2D map is formed by substituting Eq. ^ into the last term of Eq. ^ 
giving 



SXn+i 
5pn+l 



M 



6X„ 

5pn 



where 



M 



fi + vH w + vG 
H G 



(9) 
(10) 



The full control system defined by Eq. ^ has a fixed point at SXn = Spn = 0. 
The control system will produce the desired control if this fixed point is stable. The 
parameters n, v, and w are determined by the dynamics of the system near the fixed 
point and we are free to adjust H and G. For successful feedback control, the values 
of H and G must be chosen so that the fixed point of Eq. ^ is stable. 

The stability of the fixed point is determined by the magnitude of the eigen- 
values of M: 



-^1.2 



-TrAf ± V(TrM)2 - 4DetM 



(11) 



where TrM = fi + vH + G and DetM = fiG — wH. The fixed point will be stable if 



lAil < 1 and IA2I < 1. 



(12) 



Equations |Tl| and |I2| put conditions on H and G such that the control strategy will 
work. 

The best values of H and G would make Ai = A2 = and the fixed point of the 
control system becomes superstable. The superstable condition is satisfied if if = H* 
and G = G* where 



TrM 
DetM 








H + vH* + G* 
fiG* - wH*. 



The values of H* and G* that satisfy Eqs. |T3| and are 

,2 



{nv + w) 



and G* 



fiw 
{nv + w) 



(13) 
(14) 



(15) 
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Equations |T5| are identical to Eqs. ^ for the RPF algorithm with K = H* and 
R = G* and we have shown that the RPF algorithm forms a superstable control 
system. Of course, it is not necessary to have the optimum values for K and R for 
the control to be successful and there is a range of values satisfying the stability 



conditions of Eqs. O and 12. 



Finally, if we substitute H = H* and G = G* into M, then the map describing 
the superstable control system becomes 

5Xn+l = 7 r{fi6Xn + W^Spn) (16) 

[fl + W) 

= 7 7 -{-n'^SXn- fiv6pn). (17) 

(/i + W) 

Iterating these equations once shows explicitly that, for the superstable RPF control 
conditions, SXn+2 = and Spn+2 = for any starting values of {SXn, 6pn)- Of course, 
SXn and 6pn must be small enough so that the linearization of the dynamics about 
Xp is valid. 
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